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ABSTRACT 



Observations indicate that some extended outer disks have a sharp cut off in the surface 
density of neutral hydrogen when this approaches the value of ~ 2 x lO"*^^ cm~^. In this 
paper we model these HI edges as places where the ratio of neutral to ionized hydrogen 
drops rapidly due to ionizing radiation. We use two different models for the vertical 
distribution of gas above the outer galactic plane: in the first model we derive the density 
from the ideal gas law while in the second model we insert a macroscopic pressure term and 
derive the density as for an isothermal slab. We consider two different sources of ionizing 
photons: external fluxes of different intensity and spectral index due to quasars, and a 
monochromatic UV flux due to neutrino decays inside and outside the disk. We flnd that 
galaxies which have a smaller gas scale height, because of a higher dark matter density 
or a larger external pressure, should show outer HI disks to a lower column density and 
smoother HI edges. The sharpness of the HI-HII transition and the total column density at 
which the medium is 50% ionized, arc strongly correlated, irrespective of the gas or photon 
flux model used. We present several model fits to the HI sharp edge observed in the nearby 
galaxy M33. If today's UV background is dominated by attenuated quasar light, which 
gives ~ 10"-'^'* H ionizations s~^, a large gas scale height or equivalently a nearly spherical 
halo is preferred. If ionizing photons from decaying neutrinos are responsible for the M33 
sharp edge, then a thin outer HI disk, and consequently a flat dark matter halo, is required. 

Subject headings : Cosmic Background Radiation - Galaxies:ISM - Radio Sources: 21 cm 
Radiation. 



1. INTRODUCTION 



The characteristics of neutral hydrogen outside the star forming regions of disk galaxies 
have been outlined by Corbelli & Salpeter (1993), hereafter Paper I. In particular we have 
focused on the heating inputs to the gas in low pressure conditions, which favor non-LTE 
and make the HI spin temperature deviate strongly from the kinetic temperature of the gas. 
Paper I was intended for observers who have measured the HI spin temperature through 
21-cm absorption lines and want to use it to derive information on the cosmic background 
or local non-ionizing sources. In the present paper we study HI sharp edges discovered in 
spiral galaxies (M33, NGC3198) and in surveys of high velocity clouds (hereafter HVC). In 
outer disks the HI column density decreases slowly with radius until a sharp edge occurs at 
Nhi = Ni r^2x 10^^ cm-2 (Corbelh, Schneider, & Salpeter 1989; van Gorkom 1991); no 
residual HI is detected in emission at larger radii. Similarly, surveys of HVC have shown 
that there is a transitional column density AT; = 5 x 10^^ cm~^ such that appreciable 
emission is seen for Nhj > Ni but little emission is seen for lO-*^^ cm~^ < Nhj < Ni 
(Colgan, Salpeter, & Terzian 1990). In Paper I we have shown that these HI edges cannot 
be due to subthermal effects i.e. to a drastic depression of the 21-cm emission line when 
the density of the medium is lower than a critical value and the spin temperature deviates 
strongly from the gas kinetic temperature. Here we test a totally different possibility for 
HI edges i.e. that they are HI-HII transitions due to an ionizing flux which interacts with 
the neutral hydrogen distribution. 

The lack of star forming regions in the outskirts of spiral galaxies makes the cosmic back- 
ground radiation extremely important for HI ionization. This background is not directly 
observable because it is easily absorbed by the thick HI layer in the galactic plane. A 
careful comparison between sensitive HI emission studies in outer regions and theoretical 
results on the interaction between an ionizing background and a gas layer, can give im- 
portant information on the spectrum and intensity of the cosmic background below 100 
eV. Several papers have discussed this method (Felten Sz Bergeron 1969; Sunyaev 1969; 
Silk & Sunyaev 1976; Bochkarev & Sunyaev 1977); however, due to the lack of reliable 
information on low surface brightness HI and on the mass distribution in outer regions 
at that time, they could not draw any conclusion on the cosmic background. We discuss 
the properties of HI edges, how they can change from one type of galaxy to another and 
from one type of object to another (from outer disks to HVC, for example). In this paper 
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we focus on models which can reproduce the characteristics of the edge which has been 
observed with good sensitivity in M33. 

Details on the cosmic background radiation, on the equilibrium equations and parameters 
used for the gas distribution are described in Section 2. In Section 3 wc give some scaling 
relation for the gas distribution around HI edges and in Section 4 we derive numerical 
results for HI edges due to several background models and give fits to the M33 sharp edge. 
The characteristics of HI edges in outer disks which have local ionizing photons from the 
decay of dark matter are examined in Section 5. We summarize and discuss our results 
for the cosmic UV background and the gas distribution in the outskirts of disk galaxies in 
Section 6. 



2. EQUILIBRIUM MODELS FOR THE GAS IN OUTER REGIONS 

The ionization-recombination equilibrium depends on the photon spectrum and on the gas 
density distribution. The UV and soft X-ray background has been reviewed in Paper I. 
For the purpose of this paper however the spectrum above 0.5 keV is not be important 
because we are mainly interested in the ionizing effect of the cosmic radiation, and photons 
with energy close to the Lyman limit are the most important. As in Paper I we use the 
following spectral form with two parameters for the extragalactic spectrum between 13.6 
and 500 eV: 

^ = I E-^y ph cm-^s-^sr-^keV-^ (1) 

We shall mostly use a steep spectrum below 500 eV, usually s > 1.4. For s = 1.4 and 
/ ~ 7.7, eq. (1) gives the spectral law observed above 2 keV. A background dominated by 
quasar light has a spectral index s ~ 2 — 3 at zero redshift. A constraint to the extragalactic 
flux comes from Ha measurements in our Galaxy from which one can estimate the sum of 
the extragalactic and local ionizing radiation (Kutyrev & Reynolds 1989; Songalia, Bryant 
& Cowie 1989). Measurements towards HVC give the conservative upper limit of 2 x 10^ 
photons cm~^ for the total flux. The spectrum suggested by Sargent et al. (1979) for 
the extragalactic flux (7 ~ 15 and s = 2.4) gives about 5.5 x 10^ photons cm~^ and 
C ~ 10~^^ s~^ for the total number of ionizations per second for an unshielded neutral 
H atom (see Paper I, Section 3). This is much lower than the value ( ~ 2.7 xl0~^^ s~^ 
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obtained from an analysis of the proximity effect (Bajtlik et al. 1988) at redshifts of about 
2 and fits within some evolutionary studies of Lyman-a clouds (Ikeuchi & Turner 1991). 
However some recent papers on flux evolution and low redshift Lyman-a absorption lines 
estimate lower fluxes (Madau 1992; Charlton, Salpeter, & Hogan 1993). These models 
take into account the role of distant absorbers in decreasing with time the intensity of the 
cosmic radiation at energies close to the Lyman limit, and give C ~ 10~^^ s~^ at present. 
This corresponds approximately to 6 x 10~^^ erg cm~^ s~^ sr~^ Hz at the Lyman edge, 
and to eq. (1) with / ~ 1.4 and s = 2.5. In the rest of this paper we will refer to this 
spectrum as the "Madau flux" . 

The gas is considered distributed in a slab (with the vertical extension much smaller than 
the horizontal one) . For a given total column density, the penetrating flux, the gas kinetic 
temperature and Nhi /Nhii, all change with depth inside the layer. We therefore solve 
the radiative transfer, the flrst momentum and energy equations by integrating through 
the HI slab and iterate until the solution converges. Starting from the top of the slab we 
evaluate at each point the attenuation and the mean energy of the flux coming from both 
sides. We set the next step such that its optical depth, at this mean energy, is 0.05 normal 
to the plane. We use the ionization equilibrium and energy equations (eq. (8) and eq. (11) 
of Paper I) to flnd the equilibrium temperature Tk at each step. Because the mean energy 
of the radiation is close to the Lyman edge we do not include Helll and assume that the 
fractional ionization of HI {xh) is equal to that of Hel when the main photon energy is 
above 24 eV. The inclusion of Helll, as well as the correct fractional ionization of Hel, 
would make the HI edge slightly sharper because high energy photons would be absorbed 
by Hell atoms and re-emitted at lower energies with a higher probability of ionizing HI 
atoms. However, since we will be doing a best fit to the available HI data this would 
only improve the chi-square value of the fit but it will not really make a difference in the 
selection of the best-fitting model. We use half solar abundance for metals and a helium 
density which is 10% the hydrogen density. We iterate on xh and Tk at each step as well 
as on the whole solution in the slab until convergence is achieved. 

There are no observations so far of the vertical structure of the gas in outer regions (at 
radial distances greater than 2 optical radii) or in external clouds. Even if detailed studies 
of face on galaxies give evidence of broad 21-cm line profiles in outer regions (see for 
example Dickey, Hanson, & Helou 1990) we don't know yet if the linewidth is determined 
by microscopic thermal motions or by bulk motion inside the disk. Therefore we shall 
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consider separately two sets of models for the vertical gas distribution: 
Model A . This model assumes that at any height z above the gas midplane the gas pressure 
P{z) is connected to the gas mass density pi^z) simply by the ideal gas law. In this model 
we do not include the gravity due to dark matter explicitly but we allow the pressure 
above the slab, Pext = -P(oo)) to be a free parameter, somewhat similar to the parameter 
Po in Paper I. For a given total hydrogen (HI+HII) column density, Nfot-, the vertical 
distribution of the gas is then determined by the first momentum equation which in the 
absence of flow reduces to the hydrostatic equilibrium equation: 



where G is the gravitational constant. As discussed in Paper I we expect thermal pressures 
in outer regions to be much smaller than in the luminous disk and we consider 1 < Pext < 
100 cm~^ K {k is the Boltzmann constant). 

Model B . In this model we omit the external pressure, but include explicitly the gravi- 
tational force due to dark matter and allow for a contribution to total pressure P* from 
bulk motion, magnetic fields, cosmic rays, etc. We write: 



where T* must exceed the gas kinetic temperature Tk{z). While Tk{z) can vary with 
height and may change drastically from the cold to the warm HI phase, we make the 
assumption that T* can be approximated by a constant throughout the whole slab. The 
ideal gas law in this case is used only to evaluate the internal microscopic pressure. The 
density p{z) is continuos and related to the macroscopic pressure by the constant T* as in 
equation (3). Here T* is in units of (cm^ s~^) but sometimes we will express it in degrees 
K by multiplying it by k/{pmH) (with p the mean molecular weight). For the vertical 
gravity due to dark matter we use = —Kz. The exact expression for K depends on the 
halo eccentricity and we present models with different values of K. If we assume corotation 
and choose a spherical distribution of dark matter with a power law radial density profile, 
at z « R we have K ~ K^/i?^, V being the asymptotic rotational velocity. In this case 
the first momentum equation along z as: 




(2) 



P*(^) = p{z)T* 



(3) 
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In Section 1 we mentioned that, beyond the first drop associated with the optical edge, 
the neutral atomic column density decreases smoothly and slowly as one proceeds radially 
outward, until it reaches a critical value A^^ and then drops very rapidly to unobservably 
small values (Corbelli et al. 1989; van Gorkom 1991). In Paper I we showed that this sharp 
outer HI edge cannot be an artifact due to subthermal eflFects in emission, but there are still 
two different possibilities. The first one is that the HI drop corresponds to a drop in the 
total hydrogen column density. If the gas is in stable circular orbits, Ntot = Nhi + Nhii 
can be an arbitrary function and it could in principle have had a sharp drop since the 
galaxy's formation. The second possibility favored here and by Maloney (1992) is that 
Ntot continues to drop only slowly below Ni but the ionization-recombination equilibrium 
results in Nhi « Nhii for Nhi < Ni. In the two following sections we analyze the 
effects of several extragalactic background models on the ionization structure of a gas slab 
and compare this with HI observations. 



3. THE GAS VERTICAL THICKNESS AND ITS IONIZATION STRUCTURE 

For a given extragalactic fiux and total column density, Ntot, the additional input param- 
eter for Model A is the external pressure, Pext, and for Model B the constants K and 
T*. We assume that these quantities decrease smoothly and slowly with increasing radial 
distance, so that they can be considered constant around the region where the HI edge 
occurs. We then compute the vertical structure of the gas slab as we vary Ntot- We define 
the gas vertical scale height Zi/2 as the height above the midplane such that half of the 
total gas mass per unit area lies between +2:1/2 and — 2:1/2 • Similarly we define the HI scale 
height, zhi- 

For Model A the scale height depends not only on the gas surface density and on the 
external pressure Pext, but also on the extragalactic flux, since this affects Tk{z) and 
therefore the density stratification via the ideal gas law. In the optically thin regime, 
where Pext is the dominant pressure term and the temperature of the medium is uniform 
to a good approximation, the scale height is 2:1/2 = 0.5kNtotTK / Pext- For higher column 
densities, however, self-gravity is no longer negligible and the gas absorbs some of the 
extragalactic radiation. The variations of the optical depth and of the gas volume density 
along z make the kinetic temperature vary with height. The slab has a non-isothermal 
structure with a scale height which depends strongly on the gas total column density, as 
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well as on the spectrum of the incoming radiation. For this model zjji inside the HI edge is 
very small. If such thin slabs are warped or multilayered, results for edges with an external 
ionizing flux should not be changed appreciably, but observations may only give the much 
thicker envelope instead of zhi- Close to the HI edge, zhi have larger values because of 
the increased amount of warm hydrogen. 

For Model B, if dark matter dominates over self gravity, the gas density has a gaussian 
vertical distribution and ^1/2 — 0. 68 y'T*/i^ i.e. it does not depend on the column density 
Ntot but only on the ratio between the "pull down" due to dark matter and the "pull 
up" due to the inner macroscopic pressure. If instead self-gravity dominates, then the 
gas density is distributed as in an isothermal self-gravitating slab (Spitzer 1942; Ibahez 
& di Sigalotti 1984) and Zi/2 = 0.12 T* /{GmHNtot) ^ lQ^^T*/Ntot cm (where mn is 
the hydrogen mass in grams and we have used 0.1 as helium abundance in number with 
respect to hydrogen). In both cases however, the total gas for Model B has a vertical scale 
height which does not depend on the extragalactic flux but only on the local properties of 
the galaxy. The HI vertical scale height depends on the extragalactic flux but it will be a 
much smoother function than in Model A. 

For a spherical halo which has V/R ^ 6.7 km s~^ kpc~^ {K ~ 45 km^ s~^ kpc~^) the 
ratio of dark matter to gas surface density for Model B is about 4xlO^°/A^tot- The density 
distribution along z near the HI edge is close to a gaussian for most of the models we 
consider since Ntot ~ 10^*^ cm~^ and dark matter dominates over self gravity. For other 
models self-gravity is not completely negligible and we always include it in our numerical 
analysis. The observed variation of the thickness of the HI layer with radius have been 
derived for the Milky Way by Kulkarni, Blitz, & Heiles (1982) and most recently by 
Merrifleld (1992). In the outermost part {R ~ 20 kpc) Merrifleld's measure for the observed 
vertical thickness, /i^, is of order 1.3 kpc and it increases beyond as A/i^ fa O.lSAi?. If the 
HI distribution is assumed to ge Gaussian then is the rms of the Gaussian distribution 
and Zhi — O.Shz- An HI edge which occurs further out, say at about 25 kpc, will have 
hz ^ 2 kpc and zhi — 1 kpc. A spherical halo with a rotational velocity of about 220 km/s 
and 2:1/2 ~ Zhi implies for Model B a value of T* ~ 20, OOOK, or a dispersion of about 13 
km s~^. As might be expected this value for the parameter T* is slightly larger than the 
gas kinetic temperature of warm HI or than the HI velocity dispersion measured in face-on 
galaxies (Huchtmeier & Bohnenstengel 1981; Shostak & van der Kruit 1984; Dickey et al. 
1990). Since we don't know how the parameter K (which depends on halo eccentricity) 
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and T* vary from galaxy to galaxy, we shall consider models for a wide range of K/T* 
values. 

We describe the variations of Nhi with Ntot using the function: 

y = logioi^ (7) 

and define N112 as the total column density which is half neutral and half ionized, i.e. 
y{Ni/2) = 0.3. 

Consider first a hypothetical monochromatic flux of ionizing photons with photon energy 
E and ionization rate C, impinging on a layer of total column density Ntot (a variable) and 
of constant hydrogen volume density n. Let Ne be that neutral hydrogen column density 
for which the optical depth is unity, 



IE 



i^^f^)%m- (5) 



where gi/ is the gaunt factor for photoionization from ground level (Spitzer 1978). As in 
the simple "Stromgren sphere theory", we assume that the dimensionless quantity C,/{a2n) 
is large and we define a critical column density 

iVc = -^Ne > Ne (6) 
na2 

where a2 is the recombination coefficient excluding captures to the first level. The neutral 
column density Nhi then varies with Niot as follows: for Ntot > (2 or 3) Nc, most of 
the hydrogen is neutral, shielded by and outer layer with Nhii ~ Nc and y is small. For 
^tot < O.SATj,, on the other hand, the hydrogen is mostly ionized with 

AT ~ AT Ntot 

Nhi ~ = ~/V~ 

To summarize this simple model: as Ntot drops from a little above Nc to a little below, 
Nhi /Ntot drops rapidly. The larger {(/na2) is the more rapid is the drop. Observed values 
of Ni where edges occur should then be of order of N^. 

Realistic fluxes are non-monochromatic, and the shape of y as a function of Ntot depends 
on the shape of the photon spectrum. Large values of s in equation (1) should give results 
similar to that predicted by monochromatic photons at energies just above 13.6 eV; Nhi 
predicted outside the edge will be small and unobservable with present-day 21cm emission 
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techniques. On the other hand for small s (flat spectrum), the variations of y{Ntot) should 
be more gradual, since photons of different energies have different values of A^^;. 
For the variation of conditions with radial distance in a realistic disk, the density n 
varies as well as A^tot- For Model A we assume Pea;t=constant and as Ntot decreases 
radially outward, the gas pressure P in the plane decreases at first and then remains near 
Pext- The value of y increases between zero and a maximum constant value for small 
Ntot, Vmax ~ logiQ(C/nQ;2)- For Model B as Ntot decreases radially so does P* and n; y 
continues to increase slowly for small Ntot- Numerical results on the sharpness of HI-HII 
transition are discussed in the next section. 

4. NUMERICAL RESULTS USING THE EXTRAGALACTIC IONIZING FLUX 

(a) A series of models 

In Table 1 we show some numerical results using Model A (Al to A4) and Model B (Bl to 
B7). The attenuated quasar background described by Madau is used in examining cases 
with different galaxy parameters. The parameters are the external pressure Pext for Model 
A and the dark matter density plus the dispersion T* for Model B. As noted in the previous 
section, for the range of parameters of interest in Model B self-gravity hardly dominates 
over dark matter gravity. Therefore results depend mostly on the ratio K/T* (and not on 
the two parameters separately) and we show them for three values of this ratio. For each 
model we give in Table 1 the value of A''i/2, -^1/2 and zhi', for models where 2:1/2 and zhi 
vary with Ntot they are given for Ntot = 

For the models Al, A2, A3 an increase in pressure increases the mean volume density 
n, and consequently decreases the transition value A^i/2 (which is of the order of AT^, in 
eq. (6)). For a given spectrum A^^ is constant and a decrease in N1/2/NE makes the HI 
edge less sharp. If, for a given spectral shape, P and n were constant throughout the slab 
then results would depend only on the ratio I/P which would take the place of (/na2 in 
equation (6). In reality there is no exact scaling, since the internal pressure increase in 
the slab depends on Pext/A^tot (^®® only when this ratio is large the pressure is 

radially constant (P = Pext)- For the Model B series at a given spectral shape we have the 
approximate scaling A''i/2 \/Izhi and zhi ~ (2/8)2:1/2 oc a/ K/T*. For Model A zhi is 
surprisingly small for the smallest Pext because given the large value of A''i/2 a small cold 
core forms at the center of the slab. 
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We also give results for different assumptions on the background: the spectral index is the 
same as that given by Madau but the intensity is ten times stronger (B4) or the background 
has a spectral index different from s = 2.5. For model B5 we use a flatter spectrum which 
is an extrapolation down to 13.6 eV of the hard X-ray background spectrum, and for 
model B6 we use a steeper spectrum with s = 4 and a photon intensity at the Lyman 
edge ten times stronger than Madau's flux. In Figure 1 we show the total and HI column 
density as function of the dimensionless parameter x = 20— logioA^tot for various models 
listed in Table 1. Notice that Case Al and Bl give very similar HI edges despite the fact 
that they have diff'erent vertical scale heights. Futhermore the characteristics of HI edges 
corresponding to different spectra but to the same N1/2 are similar. From Figure 1 we 
note that the HI column density where the curve Nhi{x) is steepest is very small and 
unobservable with present techniques; for even smaller values of Nhi the curve Nhi{x) 
flattens out. 

(b) Fitting the HI edges observed in M33 and in HVC 

Corbelli et al. (1989) show in their flg. 10 the second HI fall-off along the north side 
of the major axis of M33. For comparison with our models, the observed values of Nhj 
should all be multiplied by an inclination factor cos(z); this factor is uncertain and we 
initially assume it is close to one. Omitting the very last point with large errors, we plot 
in Figures 2, and 3 the values of Nhi{R) at the five outermost points (filled circles). Since 
the fractional change AR/R across these points is small, we assume that Pext and K/T* 
are constant across the HI fall-off. Each model listed in Table 1 gives a unique value of 
Ntot for each observed value of Nhi- These are shown as open circles in figures 2 and 3. 
In order to compare a model with the observations we have to make some assumptions 
about the spatial variation of Ntot- We choose Ntot{R) — Nq x 10^/^' and determine 
A^o and the scale length Ri by a two parameter least squares fit for each models. We fit 
\ogiQNtot{R) assuming an RMS error of ±0.05. This is about half the error of the more 
rapidly varying logNni- The value of the fitted Ri and the x^- values for the mean squared 
deviations (with 3 degrees of freedom) are given in Table 1 for each model. 

In the outer disk of M33 the radial fall-off of Nhi changes drastically beyond 3 x 10^^ 
cm~^, suggesting that Nc ^ Ni :^ 2 x 10^^ cm~^ and N1/2 ~ 4 x 10^^ cm~^. Table 1 
and the fits in Figure 2 and 3 show indeed that models with A1/2 ~ 3 x 10-^^ cm~^ give 
a poor fit and have large values of x^- They require a sudden drop in Ntot to fit the 
sharp edge i.e. if Ntot is a smooth function of R they predict HI edges less sharp than 
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observed. These models are those with large P^xt or K/T* values (A3 and B3 in Fig. 2(c) 
and 3(c) for example) or with flat backgromid spectra (B5 in Fig. 3(e)). They are rather 
unattractive for the face-on case also because they predict very small values for Ri in M33. 
Although observations do not give a precise value of Ri, we know that Nhi{R) before the 
edge occurs drops by a factor 10 in about 20 arcmin (Ri 4 kpc, see Figure 2 of Corbelli et 
al. 1989). Therefore we don't expect Ri to be smaller than 20 arcmin for a face-on galaxy, 
improves as increses. However values of larger than those shown in Table 

1 require zhi > 5 kpc i.e. K noticeable smaller than estimated (i^(M33) > 40 km^ s~^ 
kpc~^) and unresonably high values of T* are needed (T* > 10^ km^s~^). This is to avoid 
self-gravity compression for these large column densities. 

If the outer disk of M33 is not face-on but is inclined by an angle z, the neutral column 
densities perpendicular to the galactic plane are all smaller than those observed along the 
line of sight. Assuming cosz ~0.7 or 0.5 with fluxes and galaxy parameters as given in 
Table 1, the resulting total gas scale length, i?^, becomes larger than for a face-on disk. 
The values are slightly smaller, but as the galaxy gets closer and closer to be edge-on, 
the best value of N1/2 decreases. 

With likely values of z (i ~ 45°) and a spectral shape like that for quasars or steeper (s ^ 2) 
good flts are obtained if intensity / and pressure are such that > — 2 x 10^^ cm~^. 
For the Madau fliix models with P^xt ~ 15 or T*/i^ ~ 4 kpc^ are favored since they give 
good flts to the sharp edge with resonable sets of galaxy parameters for M33. For Model B a 
value of T* ~ 100 km^s"^ (or equivalently ~ 10^ K) is a lower limit since it implies K ~ 25 
km^s~^kpc~^, too small even for a spherical dark matter halo around M33. Slightly larger 
T* and K values are preferred. These models predict reasonable scale heights, zhi ~ 1 or 

2 kpc as observed in the outermost parts of our Galaxy (Merrifleld 1992). If Merrifleld's 
observations give evidence for an envelope of an intrinsically thinner disk with warps, 
vertical oscillations or multilayered slabs (clouds), then models with zhi < 1 kpc may be 
more appropriate. 

For H VC the shape of Nhi (R) is unknown, but we know that Ni is smaller than in outer 
disks being Ni ~0.5 xlO^^ cm-^. A model between A2 and A3 with P^xt ~ 30 cm-^ K 
and Zhi ~ 0.3 kpc is then a possibility. Since HVC are located near the inner disk of our 
Galaxy, larger values of Pext (due to coronal halo gas) seem reasonable. The appropriate 
models with smaller N1/2 predict less sharp HI edges. 
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5. HI EDGES AS PREDICTED BY DARK MATTER DECAY PHOTONS 



The dark matter decay theory recently developed by Sciama (Sciama 1990, 1991), predicts 
a universe populated by heavy neutrinos which decay radiatively producing photons with 
energy = 13.6 + e with e ~ 1 eV and a production rate proportional to the density of 
dark matter. In this Section we examine the possibility that sharp HI-HII transitions in 
outer regions are due to these low energy photons generated both inside and outside the gas 
slab. If the local density of the gas is sufficiently high (as in the optical region of galaxies), 
these photons are produced and absorbed locally and they represent the most important 
source of photons for HI ionization. However, in considering the effects of neutrino decay in 
outer regions, we have to consider two other possible source terms. Due to the small value 
of Nhi, contributions to the local flux from low density regions surrounding the galaxy, 
where photons from dark matter decay are not all absorbed locally, might be significant. 
We will refer to this as the "halo contribution" . Furthermore, an additional contribution 
is given by the extragalactic UV photons created by the neutrino decay in the rest of the 
universe (here we neglect the efi^ects of an additional quasar background flux). We evaluate 
the source function, as for R/Rq = 3 in the expressions given by Sciama (1990) and we 
set e = 1 eV. In computing the interaction of the UV flux from radiative neutrino decay 
with HI at the edge of galaxies, we solve the radiative transfer equation through the slab 
taking into account the contribution of all three components to the local equilibrium. The 
external flux usually dominates in the upper layers but its exact contribution to the total 
flux in each point of the slab depends upon the external pressure and self gravity for Model 
A and on the dark matter density and velocity dispersion for Model B. 

If the medium has a scale height as in Model Bl or B2 the volume density n is quite small 
and the medium is still fully ionized for column densities Ntot ~ 10^°. This is due to 
the large number of ionizing photons or equivalently to the inequality n < ^J'^Joi^2. Until 
self gravity sets in by compressing the gas vertically there is no neutral phase. Since for 
outer disks we are interested in models in which HI is present at the level of 5 — 8 x 10^^ 
cm~^ we ignore these highly ionized disks. We restrict our attention to those cases with a 
small scale height, i.e. a larger value of K/T*, which require more flattened dark matter 
halos. This means that in order to obtain the edge at the correct column density one needs 
smaller scale heights for Sciama models than for QSO background models, or equivalently 
slightly larger values of P^xt or of KjT* . In Table 1 the last five entries show the results 
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as we vary the galaxy parameters (two for Model A and three for Model B). Figures 4 and 
5 show Nhi versus x. Models AS2 and BS3 can again be excluded because of the small 
Ri and large x^. The other three Sciama models shown in Table 1 give quite good fits 
and realistic Ri values. It is interesting that, when parameters are adjusted to give the 
same A^i/2, the predicted sharpness of the edge is similar for Model A and Model B and it 
does not differ appreciably from what we get using the Madau flux instead of the decaying 
neutrinos flux. 



6. SUMMARY AND DISCUSSION 

In this paper and in Paper I we have discussed two types of HI observations which can 
be used to constrain the spectrum of the ionizing radiation below 1 keV. Observations 
show that neutral hydrogen in outer regions of spiral galaxies is warm and in Paper I we 
have calculated the heat input required from ionizing or non-ionizing sources. We have 
shown that the ionizations from background photons, together with collisions due to self 
gravity, are sufficient to bring the spin temperature well above Tfj even if there are no 
photons below 0.5 keV. This means that Nhi does not differ appreciably from the value 
Nb, inferred from the 21-cm brightness temperature. Therefore HI edges observed in the 
brightness column density of outer disks at Nb — 2 x 10^^ cm~^ (Corbelli et al. 1989, van 
Gorkom 1991) correspond to real cut-offs in the neutral phase of the hydrogen distribution. 
The observed edge column density is larger if the outer disk is inclined. 
In order to get a sharp HI edge the intensity of the background flux needs to be stronger 
than the extrapolation at low energies of the power law observed above 2 keV. In this paper 
we have modeled the extragalactic ionizing UV flux by means of an intensity parameter 
/ and a spectral index s (eq. (1)). We have used two alternative, simplifled, models to 
determine the distributions of the gas pressure P and the volume density n as function of 
height z above the plane. They have effectively only one free parameter each. In Model 
A we assume that P and n are connected by the ideal gas law, involving the microscopic 
gas kinetic temperature T^, we do not consider dark matter gravity explicitly and use a 
compression term Pexti ^ free parameter. In this model when the gas column density 
is high enough, the cold HI phase forms near the galactic plane with much larger n and 
smaller thickness than the warm HI and HII (for this thin slab, dark matter gravity is 
indeed unimportant). In Model B we introduce a dark matter gravity term, and assume 



12 



that the density stratification is that of an isothermal slab which has some macroscopic 
temperature T* . In the range of column densities we are interested in, results depend 
mostly on the single free parameter K/T*. Model B is appropriate if outer disks are fairly 
homogeneous with weak contrast between cloud structure and diffuse medium. Magnetic 
pressure, cosmic ray pressure or turbulence allow T* to be larger than Tk- We assume that 
Pext (for Model A) and K/T* (for Model B) are constant across the radial interval where 
the HI edge occurs and we compute Nhi as a function of Ntot for different background 
spectra. 

The ranges of spectral indexes and flux densities of the ionizing background radiation 
compatible with HI edges observations are also compatible with Ha measurements and 
evolutionary studies of the background ionizing flux. The spectral shape is in principle 
important; in practice we flnd the somewhat surprising but simple result that for any 
spectral index s ^ 1.4 results depend rather weakly on the spectral shape. What is more 
important is the intensity of the flux at energies close to 13.6 eV which determines the 
ionization rate per neutral hydrogen atom, C, or the parameter defined as the value of 

Ntot where Nhi = Nhii and roughly proportional to C,/n. Spectra and models which give 
similar values of predict similar sharp edges. We have shown that Model A, Model B 
or the decaying neutrino model can all give reasonable fits to the HI edge observed in M33 
if their free parameter is such that > 2 x 10^^ cm~^. If the extragalactic flux just 

above the Lyman edge is close to what Madau predicts (s ~ 2.5, / ~ 1.4), > 2 x 10"^^ 
cm~^ implies P^xt ^ 50 cm~^ K for Model A, or zhi ~ 0.5 kpc for Model B (as for an 
almost spherical dark matter halo). For a galaxy with a modest inclination respect to our 
line of sight as N112 gets larger the fits improve. Models which give 7Vi/2 ~ 3 or 4 lO"*^^ 
cm^ fit the M33 edge considerably well and require a reasonable set of galaxy parameters. 

As Pext or K/T* decreases, the transition column density, A^i/2, and the sharpness of the 
edge increase. In the case of a quasar background, the number of ionizing photons coming 
from distant extragalactic sources is likely to be constant. Galaxies with a higher pressure, 
i.e with a substantial amount of coronal gas or dark matter in outer regions, should then 
present a smaller value of A^i/2 and smoother edges. Similarly, HI column densities smaller 
than A^i/2 can survive in HVC if they have a vertical scale height smaller than that of outer 
disks due to stronger compression (from P^xt or local dark matter). With the observational 
sensitivity available at the moment, this may have led to a selection bias: edges which occur 
at larger HI column densities are easier to observe and the observed sharpness may then 
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be larger than the average. The galaxy NGC 3198 shows an HI edge rather similar to that 
of M33. This is not too surprising since they might have similar dark matter densities (for 
spherical halos V/ R ^ 5 — Q km/s kpc~^ for both galaxies). However if the dark matter 
density changes from one type of galaxy to another, then we expect corresponding changes 
in the HI edge. 

If decaying neutrinos provide the UV ionizing photons then the HI scale height must be 
quite small {zhi < 0.5 kpc) for the gas to stay mostly neutral at Ntot ~ 10^° cm~^. 
However once this small scale height is provided, by using for example a flat dark matter 
halo or a low T*, the sharpness of the HI edge is rather similar to that predicted by a 
quasar background flux which gives the same value. For the outer disk of our Galaxy 
Zhi > 0.5 kpc, however there are no observations of the HI vertical scale height just inside 
a sharp HI edge. These observations would be of great interest for selecting the most 
appropriate model (decaying neutrinos. Model A or Model B) for outer disks with a sharp 
HI edge. We hope that in the near future detailed observations of other HI edges will 
become available together with measurements of the HI spin temperature and its vertical 
scale height just before the edge occurs. 
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FIGURE CAPTIONS 



Figure 1. Ntot (solid line) and Nhi (dashed lines) as functions of the dimensionless 
parameter = 20 — logi^Ntot for models as given in Table 1. The filled triangles 
indicate the HI column density at Ntot = -^1/2- In ^(fl) we show cases for Model A 
and in 1(6) for Model B. The small dashed line refer to model Al in 1(a) and Bl in 
1(6), the medium dashed to model A2 and B2, and the large dashed to model A3 and 
B3. The dot-dashed line refer to model A4 and B4 in 1(a) and 1(6) respectively, and 
have a different background intensity than used in A1,A2,A3. In 1(c) the spectral 
slope differs from 2.5 and the models shown are B5 (large dashed line), B6 (medium 
dashed line) and B7 (dot-dashed line). B6 and B7 have the same 

Figure 2. The values of Nhi observed near the edge in M33 (filled circles) and the 
corresponding total gas density Nt^t (open circles) for a given Model A. Models 
shown are labeled as in Table 1. The straight line is the best linear fit to \ogNtot{R)- 

Figure 3. The values of Nhi observed near the edge in M33 (filled circles) and the 
corresponding total gas density Ntot (open circles) for a given Model B. Six models 
are shown and labeled as in Table 1. The straight line is the best linear fit to 
\ogNtot{R). 

Figure 4. As in Figure 1 we plot Ntot (solid line) and Nhi (dashed lines) as functions of 
the dimensionless x. Here we used a flux as given by the neutrino decay at 3 solar 
radii (internal and external). The filled triangles indicate the HI column density at 
Ntot = -^1/2- In (ti) we plot the HI column density for Case ASl (small dashed line) 
and AS2 (medium dashed line) as given in Table 1. In (6) we plot the HI column 
density for Case BSl (small dashed line), BS2 (medium dashed line) and BS3 (large 
dashed line). 

Figure 5. The values of Nhi observed near the edge in M33 (filled circles) and the 
corresponding total gas density Ntot (open circles) for four decaying neutrinos models 
described in Table 1. 
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